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ABSTRACT 

In this paper we consider Roche accretion in an Extreme Mass-Ratio Inspiral 
(EMRI) binary system formed by a star orbiting a massive black hole. The ultimate 
goal is to detect the mass and spin of the black hole and provide a test of general 
relativity in the strong-field regime from the resultant quasi-periodic signals. Before 
accretion starts, the stellar orbit is presumed to be circular and equatorial, and shrinks 
due to gravitational radiation. New fitting formulae are presented for the inspiral 
time and the radiation-reaction torque in the relativistic regime. If the inspiralling 
star fills its Roche lobe outside the Innermost Stable Circular Orbit (ISCO) of the 
hole, gas will flow through the inner Lagrange point (LI) to the hole. We give new 
relativistic interpolation formulae for the volume enclosed by the Roche lobe. If this 
mass-transfer happens on a time scale faster than the thermal time scale but slower 
than the dynamical time scale, the star will evolve adiabatically, and, in most cases, will 
recede from the hole while filling its Roche lobe. We calculate how the stellar orbital 
period and mass-transfer rate will change through the "Roche evolution" for various 
types of stars in the relativistic regime. We envisage that the mass stream eventually 
hits the accretion disc, where it forms a hot spot orbiting the hole and may ultimately 
modulate the luminosity with the stellar orbital frequency. The observability of such a 
modulation is discussed along with a possible interpretation of an intermittent 1 hour 
period in the X-ray emission of RE J1034-I-396. 

Key words: X-rays: binaries, stars: winds, outflows, accretion, radiative transfer, 
galaxies: active, galaxies: individual: RE J1034+396 



1 INTRODUCTION 

The measurement of black hole parameters - mass and spin - 
has been an important astrophysics research topics over the 
past decade. Such measurements can lead to tests of general 
relativity in the strong-field regime, and reveal information 
on the formation and evolution of galaxies. There has been 
great progress on this topic using various approaches. The 
mass of a black hole is more easily measured since it is a 
Newtonian parameter. However, the spin, a relativistic pa- 
rameter, is much harder to be measure. One method is to 
use radiation from the accretion disc: either by fitting the 
shifts and shapes of the iron fines (e.g.. iRevnolds &: NowakI 
20031) V or by usi n g a multi-temperature blackbody model 



Yuan et al.l |2004| : iLi et al.l l2009l : IShcherbakov et al.1 12010| : 



IGou et al.ll20Tl 



The problem with the iron line model is 
that it is not certain that the Fe lines indeed come from near 
the Innermost Stable Circular Orbit (ISCO). The continuum 
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method can be employed to obtain the spin of stellar mass 
black holes effectively. However, for Super-Massive Black 
Holes (SMBHs), the radiation from the inner disc region 
is in the UV range and can be absorbe d. Jets are commonly 
associat ed with black hole s p ins ( e.g., Bla ndford &: Znaiekl 
(197';^ : iTchekhovskov et al.1 l(2oTl. )). However these mod- 
els are more complicated. Here we propose another way of 
detecting SMBH spins using tidal effects induced by star- 
SMBH close passages. 

There are two extreme cases of such tidal effects. The 
first case is a tidal disruption. This can happen when a star 
passes close by a black hole in an unbound orbit or a very 
eccentric orbit. The stellar volume exceeds the Roche vol- 
ume at periastron and th e star is torn apart and eventu - 
ally swallowed by the hole (jReedl 19881 : [stone fc Loebll201ll ). 
One possible instance of a tidal disruption i s the recent Sw 
J1644 4-I-57 observation (|Bloom et al.l l201ll : IZauderer et al.1 
These events are very luminous but hard to model. In 
the other star orbiting the hole in an equatorial cir- 

cular orbit spirals inward under the action of gravitational 
radiation. When the star fills it Roche Lobe, mass is stripped 



2 Lixin ( Jane ) Dai, & Roger Blandford 



stably from the star on a time scale of ~ 10, 000 — 100, 000 
years. This steady mass-transfer is not as luminous as tidal 
disruption, but lasts for much longer. Also the second case is 
easier to model with high precision for the mass and spin of 
the hole. In this paper we will provide techniques for doing 
so. 

Mass-transfer in binary systems has been well studied 
for a variety of stellar components and mass ratios. This in- 
vestigation is mostly focused on binary systems composed 
of a stellar-mass co mpact object a nd a companion star, 



g. IWeb bink' (l9 8^). In particul ar, iHiellming fc Webbinlj 



and Sobcr man et al.l l|l997l ) explored the stability of 
adiabatic mass-transfers based on polytropic stellar mod- 
els, and showed it could be dynamically unstable. However, 
these studies are Newtonian calculations, since the compan- 
ion star reaches the Roche limit far away from the ISCO 
of the stellar-mass compact object. In Extreme Mass-Ratio 
Inspiral (EMRI) systems, stars can get close to the ISCO 
before being tidally disrupted. Therefore, the signals from 
such binary systems can in principle be used to constrain 
the mass and spin of SMBHs, and to test general relativity 
in the strong-field regime. In this paper, we present an accu- 
rate general relativistic framework for studying such mass- 
transfer. 

We shall be especially concerned when the mass-transfer 
is "adiabatic" - in other words it happens on a time scale 
short compared with the thermal (Kevin-Helmholtz) time 
for the star yet long compared with the dyna mical time. This 

f roble m has previously been considered by iHameurv et al] 
19941 ) in connection with the (erroneous) iden tification of 
periodicity in NGC 6814 (iMadeiski et all 1 19931 ). This pa- 



per builds upon the iHameurv et~ (1994) discussion and 
presents a more comprehensive treatment of stars and plan- 
ets. In particular, we will show for various stellar types how 
the stellar Roche volume and evolution change with the spin 
parameter of the hole and also the rotational angular veloc- 
ity of the star in the relativistic regime. We will also con- 
clude that such accretion events are stable o ver a dynamical 
time s cale, in contrast to the suggestions of IHameurv et al] 
l|l994h . The differences in results are mainly caused by more 
accurate calculation of the relativistic Roche potential (and 
the resultant Roche volume), as well as more accurate grav- 
itational radiation power calculation. 

In an EMRI system, the central black hole is so massive 
that its mass is effectively constant. The stellar orbit, on the 
other hand, shrinks due to the loss of energy and angular mo- 
mentum through gravitational radiation and other torques. 
This radiation is a prime target for space-borne gravitational 
radiation observatories intended to test the general theory 
of relativity in the strong field regime. An EMRI system can 
also emit electro-magnetic radiation (principally X-rays and 
optical), and its observation can also provide the basis of a 
good test of general relativity, if measurements are carried 
out accurately. 

We restrict our discussion of EMRI systems to stars in 
circular, equatorial orbits. If and when such a star overfiows 
its Roche lobe close to the ISCO, gas will flow through the 
inner Lagrange point (LI), and the stellar orbit will usu- 
ally expand - an Extreme Mass Ratio Outspiral or EMRO! 
We shall assume initially that the stellar evolution is adia- 
batic, and discuss later the limitations of this assumption. 
Under these conditions, the entropy variation with mass 



within the star is frozen and the evolution is as described 
in our previous paper. In the Newtonian regime it is the 
sign of {dp-i,/dm^)s, where is the total remaining mass 
of the star, and p is the mean stellar density, that determines 
whether or not an EMRI becomes an EMRO. 

The gas flowing through LI will generally have enough 
angular momentum to form a torus orbiting the hole. A hot 
spot will be formed where the stream hits the disc. In ad- 
dition, it is quite possible that one or two spiral arm shock 
waves will be created in the disc. The hot spot and shocks 
will orbit with the orbital frequency of the star and an in- 
clined observer should see a modulation of the radiation 
from the disc, most likely in the X-ray range. Clearly the 
modulation will be the greatest when the observer is in the 
equatorial plane, observing maximal Doppler boosting. 

One observation that may possibly be explained by this 
model is the quasi-periodic XMM-Newton X-ray signal from 
a nearby active S eyfert I galaxy RE J 1034-1-396 observed by 
ICierlinski et al" I (|2008t ). As shown in Fig. 1 in their paper, 
the rest-frame X-ray emission has a periodicity of ~ 1 hr. 
The X-ray flux has a modulation of ~ 0.1 associated with 
the nucleus of the galaxy, and this modulation could be well 
fitted by a cosine function. We will show in this paper that a 
dwarf star or brown dwarf or red giant can produce emissions 
with such a periodicity filling a ~ lO'* - lO^M© SMBH un- 
der Roche accretion, though the mass accretion rate under 
adiabatic condition and pure gravitational radiation torque 
is too low to produce the observed X-ray flux. 

In section 2, we will illustrate the mechanism by com- 
puting the evolution in the Newtonian limit. In particular, 
we will show how the star spirals in, when it reaches the 
Roche limit, and we will calculate how the stellar orbital 
period and mass-transfer rate evolve taking a dwarf star as 
representative. We will also discuss the final fate of stars 
after they have undergone Roche evolution. In section 3, 
these calculations will be repeated in the rela tivistic regime, 
and be compared with lHameurv et al.l (|l994l )'s. In section 4, 
we will use the data from RE J 1034-1-396 to fit our model, 
and to make predictions of the future signals. In particular, 
we will discuss how spins can be constrained from signals. 
Discussion and future research directions will be laid out in 
the last section. In our forthcoming paper in preparation, 
we will use these new results to construct new models of 
EMRI/EMRO evolution. 



2 NEWTONIAN TREATMENT 

A Newtonian trea tment brings out the salient features of 
the problem, e.g., IHameurv et al.l (|l994l ). Consider a fully 
convective lower main sequence star with a mass of m* = 
mMo as a representative, and assume that it orbits a hole 
with a mass in a circular, equatorial orbit 

with radius r. (The stellar orbit circularizes quickly either 
through gas dynamical drag or gravitational radiation as 
shown in iDai fc BlandfordI l|201ll ). We will take G = c = 1 
henceforth. 

From Kepler's third law, the orbital radius (measured 
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in units of the gravitational radius Rg{Ale) = 1.47 x 
10" Me cm) is: 



f = — = 23.9 P3Mg ^ (1) 

Hg 

where the orbital period P is measured in hours. 
The orbital angular momentum: 

L = m*VMr = 4.9 mM^ pi (2) 

measured in units of Mq c Rg(Ma = 1) = 8.8 x 10^'*g cm^/s. 

Likewise, the orbital energy, in units of Mqc? — 1.79 x 
10^* erg, is 



E = -0.5A/m*r"' = -0.021 x m.M§ P~^ . 



2.1 Pre-Roche evolution 



(3) 



We assume that the star follows a circular orbit which 
shrinks slowly due to gravitational radiation. In the New- 
tonian limit, the torque is 

32,«-| 2/27r 7 44 2 ,^1 r.- i 2-2 

T = —M^m* (-5-) -^ = 1.7 X 10 m Mq^ P cm s . 

(4) 

The time it takes the star starting from radius r to 
plunge into the hole is (assuming this gravitational radiation 
torque, and r ^ nsco, which is the radius of the ISCO, see 
below): 



^inspiral — 



-2 s 



' 987 m"'Mg yr. 



(5) 



256 M^m* 

(e.g. iMisner et all (|l973l )). 

Typically if the orbit is at least mildly relativistic, the 
inspiral time is less than a Kelvin-Helmholtz time. 

2.2 Roche limit 

In order to obtain the Roche volume of the star, it is con- 
venient to transform from an inertial frame to a frame of 
reference rotating with the binary system. 

The Roche potential <I>(r) for an arbitrary point with 
position vector r in the rotating frame is given by: 



near the star, and note that the linear term is canceled 
by a linear term in the expansion of the centrifugal po- 
tential, the Newtonian tidal tensor contributes quadratic 
terms —^{cf>,xxx'^ + (t>,yyy'^ + (t>,zzZ^), where (t>,xx = —'2(t>,yy = 
—20.22 ~ 2M/r^. (Note that the contribution to the Lapla- 
cian vanishes as it should.) The expansion of the centrifugal 
term also contributes a quadratic form ~^lj^{x'^ + z^), or 

A co-rotating star will expand to fill an equipotential 
surface. However, it is rather unlikely that an inspiralling 
star will be able to maintain corotation. If it rotates with 
angular frequency ui < Q,, the quadratic centrifugal term 
is approximately —^u'^ix'^ + z^). which includes the effects 
of Coriolis force. (There may be some internal circulation 
which will vitiate adiabaticity, but we do not consider this 
possibility further in this paper.) 

The gradient of the Roche potential vanishes at the first 
and second Lagrange point. Using this we can calculate for 
two extreme cases that the first and second Lagrange points 
(LI and L2) are located at 



XR = ±0.794 X 



(7) 



from the star, where X{kuj) decreases from 1 to 0.87 as fc„ 
goes from to 1. 

The transverse extents of the equipotential Roche sur- 
face are: 



y« = ±0.497 3^(fcj(^)- 
and 



ZR ^ ±o.m z{k^) r 



(8) 



(9) 



where y{kui) goes from 1 to 0.89 and Z{kui) goes from 1 to 
0.93 as kui goes from to 1. 

Furthermore, the volume contained by the Roche sur- 
face is numerically calculated to be: 



... 0.683 x(^±l)" 



(10) 



Combining this with Eq. [T], we obtain the relationship 
between the orbital period and the mean density of the star 
when the star reaches the Roche limit: 



$Miv = -— — -— ^---c3xf^' 6 
|r--r*| |r-rBH| 2 

iFrank et all (|2002| ). where r* and rsH are the position vec- 
tors of the centers of the star and the hole, and uj is the 
rotational angular velocity of the star in the inertial frame. 
uj satisfies Q ^ uj = k^Q, < Q,, where Q. is the stellar orbiting 
angular velocity. 

The first term in this equation in the Roche potential is 
the gravitational potential of the star itself, making the ap- 
proximation that the stellar mass is centrally concentrated 
and only retaining the monopolar term. The second term is 
the gravitational potential term of the hole. The third term 
is the potential in the rotating frame due to the centrifugal 
force. 

We erect a locally orthonormal coordinate system 
{x, y, z} centered at the star and directed along the {r, 6, (j)} 
coordinate axes. If we expand the gravitational potential 



Po^s.ml^^ + iy -p-'/' (11) 

independent of the mass of the hole, where Pq is the stellar 
orbital period in hours, and p the mean density of the star in 
cgs unit. For a dwarf star with mass ~ 0.3Mq , the Roche 
condition is satisfied when the orbital period decreases to 
a value of 1.96 — 2.27 hr depending on k^j. If we suppose 
that the star plunges when its orbital radius reaches 6Rg, 
the maximum black hole mass for there to be mass-transfer 
from a ~ 0.3Mq dwarf star is 1.6- 1.9 x lO^'M© {k^ = 0-1). 

2.3 Roche evolution 

To set the stage for the relativistic treatment below, let us 
consider the classical Roche model for a dwarf star that is co- 
rotating. Its mean density decreases monotonically as mass 
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is tidally stripped, leading to expanding orbits. Assuming 
that the star has initial mass m^o ~ O.3Af0, we found that 



(12) 



ijPai fc Blandfordl(|201ll )). where p*o is the initial mean den- 
sity of the star. 

Combining Eq. [11] and [12] , the orbital period P of the 
dwarf star, in terms of the initial period Pq, varies with its 
remaining mass as: 

(13) 

Po \m^o J 

We assume that the orbit of the star will change in a 
way that the angular momentum of the system is conserved 
so that 

.+ f^0, (14) 

where we continue to assume that gravitational radiation 
provides the only torque as in Eq. [1]. 

Using Eq. [2] , [4] , [13] , and [14j , we see the time elapsed 
since the onset of mass-transfer goes as: 

3.6 



tR ~ 4.1 X m-'M^ ^Po 



-V 



1 yr. 



(15) 



The associated rate of mass transfer of Roche accretion, 
also obtained assuming conservation of angular momentum, 
is: 



-1.23 X 10^^ A'/, 



6 



-3.62 X 10^'' Po 



-2. 7,, 0.7 4.9 -1 

Me m g s . 



(16) 



We show in Figs. [T] and [2] the Newtonian evolution of 
this low-mass star in orbit around a hole with mass IO'^Mq. 
We see that the stellar orbital period can increase to around 
five times the Roche period as the star is stripped. How- 
ever, the majority of the Roche evolution happens on a time 
scale ~ 10^ years, and the stellar mass-transfer rate deceases 
rapidly with time. 

As discussed by iDai fc BlandfordI (|201ll '). the angu- 
lar momentum of the binary system as in Eq. [2] follows 
1 1 

L ~ m-^Ps ~ mi,Pi, " . If the mass accretion is stable on 
dynamic time scales, it mean L should decrease as mass 
is accreted onto the hole. If we define a parameter rj = 
d(lnp*)/d(lnm*), then ?7 < 6 is the criteria for Newtonian 
dynamical stability. 

A more intuitive method of c alculating the dynamical 
stability, as .Hameurv et al.l (|l994l 'l suggested, is to ensure 
that the star will not exceed the Roche surface under the 
assumption of no angular momentum loss through gravi- 
tational radiation. In other words, dv^/dm^ > dvt^/drrii,, 
where = rri-^lpi, is the stellar volume, and un is the Roche 
volume in Eq. [10]. This leads to the same stability criteria 
r; < 6. 

Continuing to consider a ~ O.3Af0 dwarf star, we have 
r\ — 2.2 (Eq. [12]), ensuring that the star recedes from the 
hole stably under Roche accretion. 

For many purposes this Newtonian treatment applied 
to a broader range of stars will be adequate. However, as 
we want to exhibit relativistic effects, we will turn to a rel- 
ativistic treatment, a gain building upon the discussion in 
iHameurv et"al] (|l994l l. 
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Figure 1. This figure shows the period-mass relationship of a 
non-rotating lower main-sequence star with initial mass ~ 0.3Mq 
accreting onto a Schwarzschild SMBH with a mass of IO'^Mq. The 
vertical dashed line represents the evolution prior to the mass- 
transfer, when the mass does not change, but the orbit shrinks. 
The black curve describes the subsequent Roche evolution in the 
Newtonian limit, when the stellar orbit expands. The blue curve 
shows how stellar orbital period changes with mass in the rela- 
tivistic limit. 
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Figure 2. This figure shows how fast a non-rotating lower 
main-sequence star with initial mass ~ O.3Af0 accretes onto a 
Schwarzschild SMBH with a mass 10^ Mq through the Roche evo- 
lution. The black and blue curves are Newtonian and relativistic 
calculations correspondingly. 



2.4 The final fate of stars under Roche evolution 

After mass-transfer begins, most stars will recede from the 
black hole so that they will continue to fill their Roche 
lobes as their density falls. The star loses mass less and 
less rapidly. Eventually the star will cool and its structure 
will evolve. The motion of the star under the action of grav- 
itational radiation will then cause the orbit to shrink. Ad- 
ditional excursions are possible depending upon the struc- 
ture of the star. The star might cool down adiabatically or 
on a thermal time scale. Therefore, it might or might not 
fill it Roche Lobe on its return. The equilibrium proper- 
ties of these stars will be even more complicated than the 
ones we calculated in this paper. Ultimately the star crosses 
the ISCO and plunges into the hole. When more than one 
star is in a "parking orbit" there is the possibility of reso- 
nant dynamical interaction analogous to that which occurs 
in planetary systems. In addition interaction with the disc 
can alter the evolution. 

We illustrate the scenario in the following Fig. |3] Around 
the hole, there will be stars migrating inward and mostly 
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Figure 3. For a ~ Mq star orbiting a Sciiwarzschild SMBH witii 
mass IO^Mq as discussed, tlie star readies tlie Roche limit when 
the orbital radius is ^ 12Rg, and then spiral out. When the evo- 
lution time scale is comparable with Kevin-Helmholtz time scale, 
the star stops receding in the "parking belt", which is located 
at around ~ 25-Rg from the hole. Here we calculated the Kevin- 
Helmholtz time scale of the star using the initial stellar internal 
energy and luminosity. 



outward under Roche evolution, stars that have gone at least 
one Roche evolution "parked" in a belt or cloud (if the stars 
came in on inclined orbits) waiting to move in again, new 
stars coming in under gravitational radiation still burning, 
and stars coming in from the "parking cloud" and not burn- 
ing. 



3 RELATIVISTIC TREATMENT 

3.1 Gravitational radiation emission 

A star (with a mass m*) orbiting a massive hole (with a mass 
M and a spin parameter a) will lose angular momentum and 
energy due to gravitational radiat ion. The orbit is c i rcular - 
ized quickly, as we have shown in lDai fc BlandfordI l|201lD . 
However the rate of alignment of the spin and orbit al angu- 
lar mo mentum is slower than the inspiral rate fe.g.. iHughed 
For the moment we restrict attention to circular 
equatorial orbits. 

We define a dimensionless orbital angular velocity f2 = 
MQ, where Q is the regular orbital angular velocity. The 
relativistic Kepler's law is: 



f 2 rt a 



(17) 



where ± refers to prograde and retrograde orbits. 

The star plunges into the hole when the orbit shrinks 
to the ISCO. The ISCO of a black hole of mass M and spin 
parameter a has a radius risco that satisfies the equation 
as given m iBardeen et all (| 19721 ): 



.CO 0.0 
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Arctan(r//fisco) 



Figure 4. The relativistic correction term T to the plunge time 
as a function of the distance to the hole and the spin parameter 
of the hole. We see that as the radius goes to infinity, T becomes 
1 no matter what the spin is. Also 7" = when r = -Risco- 



This radius equals Rg, 6Rg, and 9Rg for a — 1, 



and — 1 

respectiv ely. The upp er limit for a has been argued to be 
~ 0.998 (|Thornelll973 ) for which nsco ~ 1.237Rg. 

Let us now consid er how to modify the treatment of 
gravitational radiation. iFinn fc Thornd l|2000l ) computed a 
relativistic inspiral time scale for a star at an orbital radius 



ti 



nspiral 



256 mMi 



T 



(19) 



in the Boyer-Lindquist coordinate, where T is the relativistic 
correction term. They presented their results of T in tabu- 
lar form. We have used these results to give a convenient 
interpolation formula for T: 



T{x, a) 

b 
d 
f 



{x- if + b {x - 1) 
x'^ + d X + f ' 

0.11 



(20) 



1 - ay 



1 - a)"'^^ - 2.22, 



1.06 X 
0.57 X 
hl.OSx (l-a)^-""- 

^g-0.23x{l-a)^4^g. 



0.66, 



3.18 X (1 - ay 



where x = r/risco- Clearly T{x — l,a) = 0, and T{x,a) — >■ 
1 as a: — ^ -l-oo (since at infinity the relativistic effects are 
negligible). This formula holds true with a maximum error 
of ~ 2% for a; ^ 2. A contour plot of T{x, a) is shown in 
Fig.H 

The results for the evolution of the 0.3Mq dwarf star 
are shown in Fig. [S] It is clear that the standard "Landau- 
Lifschitz" form equation is adequate for r > 100-Rg . 

The gravitational radiation reaction torque is readily 
computed from: 



' nsco 
. M , 



, nsco 



-8a 



nsco „ 2 

-nrr 

M 



(18) 



dtpi 

un 

dr 



(21) 
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Figure 5. The total inspiral time in years as a function of r/Rg 
for a 0.3Mq dwarf star orbiting a IO^Mq iiole. Tliis plot is in 
the log-log scale. The three colors indicate three different spin 
parameters of the hole: -0.99 (green), (black), and 0.999 (blue). 
The three solid curves are accurate inspiral times obtained using 
Eq. 1191 . The red dotted curve is the Newtonian inspiral time 
using Eq. [S]. All other curves approach the Newtonian curve 
when the star is far away from the hole. 



where the orbital angular momentum 

L = m^M^ (22) 

(1 - 3f-i -I- 2a* f"2)2 

fe.g- lBardeen et al.1 l|l972l )). 
3.2 Roche Hmit 

We use the same approach as the Newtonian treatment. The 
local tidal tensor is evaluated in Boyer-Lindquist coordinates 
and then transformed into an orthonormal coordinate basis 
freely falling with the star. This "local" coordinate system 
has spatial coordinates directed along the {r, 6, cj>} axes. The 
only non-zero spatial elements of the tidal potential can be 
expressed in the form: 

4>,xx = (P,yy = — , 0,zz = — , (23) 



where 



k = 



4af 2 -I- Sa-^ 



(24) 
— 3f + 2af 2 

This satisfies Laplace's equation as it must. Note the re- 
markable prospects that fe — > 2 for r — >■ risco independent 
of a. Also fc — >■ 1 as r — >■ 00 as must also be true in the 
Newtonian limit. 

Combining the stellar gravity and rotation with the 
tidal contribution, the Roche potential becomes: 



/ 2 , 2 

— m*(a; +y 



+ 



M ky^ + -{k + l)x^ 



Mk^{x^ + y^) 




Figure 6. Contour plot of Roche volume as a function of k and 
fcu for a = 0.5. A; = 1 is the Newtonian limit, k = 2 when the star 
is at ISCO. k^ = when the star is non rotating, and k^ = I 
when the star is co-rotating. 



Here k^^ is the stellar rotational parameter in the inertial 
frame. 

Similarly to the Newtonian treatment, we can locate the 
first and second Lagrange points and the transverse extent 
of the Roche surface, and calculate the volume contained by 
this surface. When k — 2, the Roche volume at the ISCO is 
calculated to be: 

-1 



Visco = 0.456 X 



4.09 



+ 1 ^ 



m^, 3 



(26) 



And more generally, an adequate approximation to the 
Roche volume is: 



V{r,a,ku;) = VNir,k^)- 



V^ir, k^) - Visco(r, k^ 



(^— ) +F(a,fe„)(- 



-1) 
(27) 



where 
F(a, k^ 



-23.3 + 
+ (^23.8 
+ ( 0.9- 



13.9 

k^ -f 2.8 
14.8 



2.8 + K 
0.4 
2.6 -f ku, 



X (1 - a)~ 



(25) 



As an example. Fig. [5] is the contour plot of the rela- 
tivistic Roche volume (in units of m*/Mr'^) around a black 
hole with a = 0.5 as a function of k and k^j,. This volume 
decreases from 0.683 to 0.456 from the Newtonian limit to 
the most extreme relativistic limit (r = nsco) for a non- 
rotating star. For a co-rotating star, the variation is from 
0.503 (fc = 1) to 0.368 (fc = 2). These volumes at the four 
corners of the plot ar e independent o f the spin. 

Comparing wi th Fishbonj l| 19731 ) 's Roche volume used 
bv lHameurv et al.l l|l994^ ■ our model is more accurate in the 
relativistic limit, including black hole's spin and the stellar 
rotation. We plot the ratio of two Roche volumes in Fig[7] 
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1.4 




0.8 



20 40 60 80 100 

r/R, 

Figure 7. The ratio of our Roclie volume to iFishbond l|l973l )'s 
Roche volume, as a function of r/Rg. The our four solid curves are 
for non-rotating stars orbiting a black hole with a = —I (red), 
(blue), 0.9 (black), and 0.998 (green). The lower four dashed ones 
are for co-rotating stars correspondingly. 

Using Eq. [T7] and Eq. [57], we obtain a relationship 
between the stellar mean density and the orbital radius when 
the star fills its Roche lobe: 

M _ 0.683 X -fl) ~' - 0.456 x ( + l) ~' 

(-^Y'" + F{a,k^){—^ 1) 

V nsco / V ' / \ nsco ' 

+0.683 X + \ (28) 

Then we know how the stellar properties govern the orbital 
evolution during the Roche accretion phase. 

3.3 Roche evolution 

We studied h ow a star evolved adiabatically in 

iDai fc BlandfordI ^201^ . and calculated how the stel- 
lar mean density changes as mass is stripped for a 
representative set of stars. As the star approaches a SMBH 
in a circular, equatorial orbit, it might plunge into the 
SMBH directly, go through full Roche evolution, or start 
the Roche evolution then plunge into the hole. The detailed 
evolution depends on the mass and spin of the black hole, 
the stellar type, and the rotational parameter fc„ of the star. 
Using Eq. [28] , we can check when the star will start Roche 
accretion and if so, how its orbit will change through the 
process. We summarize these inspiral/outspiral scenarios in 

Each plot shows how different stars behave as they spi- 
ral to within 50-Rg of a SMBH with given mass and spin. 
General relativity plays an important role in this region. 
For example, in plot (f), the non-spinning SMBH has mass 
IO^Mq. Brown dwarfs and white dwarfs will plunge into the 
hole. A Sun-like star, an upper main-sequence star with an 
initial mass ~ 7.9Mq, and a lower main-sequence star with 
an initial mass ~ O.3M0 will go through full Roche accre- 
tion. A red supergiant with an initial mass ~ ll.SSM© will 
start Roche accretion and will plunge into the hole. 

The plots also show how the stellar orbital periods 
change in the Roche phase. The Period-Mass curves can 
be quite distinct with or without relativity or if the star 
is CO- rotating with the black hole or not. In plot (c), a Sun- 
like co-rotating star barely crosses the ISCO of a W^Mq 



M = IO'Mq , a = 0.9, P,s<.„ = 2.04 hr 



50 



















Jupiter \ \\ 










Earth 





10-' lO-' 10-^ lO-" 0.001 



Figure 8. Roche evolutions of a perfect solid planet with the 
mass of Earth (green) and a perfect gas planet with the mass 
of Jupiter (pink) near a maximally spinning SMBH with M = 
10* M0. We plotted four curves for each planets for comparison: 
the relativistic calculation for a co-rotating star (colored solid 
curve), the relativistic calculation for a non-rotating star (colored 
dashed curve), the Newtonian calculation for a co-rotating star 
(corresponding black solid curve) , and the Newtonian calculation 
for a non-rotating star (corresponding black dashed curve). The 
vertical blue dashed lines represent the phase of inspiral due to 
gravitational radiation. 

SMBH with a — 0.9 using Newtonian calculation so it can 
go through Roche evolution. However, using relativistic cal- 
culation this star will just plunge into the hole. And in plot 
(1), only with general relativity a co-rotating white dwarf 
with initial mass O.8M0 can start the Roche accretion phase 
near a maximally spinning M6 black hole. 

We also plot the Period-Mass relationship for planets as 
in Fig. [8] A solid planet like Earth has the same period as 
mass decreases due to unchanging density. A gas giant like 
Jupiter evolves so that the period will increase. 

From the conservation of angular momentum, the stel- 
lar angular momentum L satisfies: 

AL = ALgr + ALisco + AL„ind + AL.ad- (29) 

Here L is the angular momentum of the star as in Eq. 
|22) . ALgr is caused by the gravitational radiation torque 
Tgi- calculated from out new model (Eq .[21]). ALisco repre- 
sents the loss of angular momentum of materials falling into 
the ISCO. ALwind is the angular momentum of materials 
lost from an accretion disc. In this inner region of an AGN 
accretion discs, there are mainly two kinds of wind: radia- 
tion pressure and magnetic driven. The radiation pressure 
induced wind carries its own specific angular momentum as 
on the disc. The magnetic induced wind carries away very lit- 
tle mass but very large angular momentum which can be 10 
times as large as its specific angular momentum on the disc. 
The last term is the disc radiation torque due to beaming 
effect, which is relatively small compared with the others. 
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(a) M = 10"A/q, a = -1, = 22.35 hi 



(bl M = ID'A/g , a = 0, Pi„„ = 12.64 hr 



(t) M = lO^A/g , a = 0.9, Pi^c = 3.81 hi 



(d) M= IO^A/q. a = 0.998, Pi„„ = 2.04 hi 
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Figure 9. Roche evolutions of different stars close to SMBHs with M = 10* M©, IO'^Mq, lO'^M©, IO^Mq, and a = -1, 0, 0.9 and 0.998. 
In each plot, the x-axis is the stellar mass ranging from O.OOSM© to ISM©, and the y-axis is the stellar orbital period from Pisco to 
the orbital period at 50-Rg (close enough to the SMBH so these are the relativistic regions). We use different colors to denote different 
types of stars: a Sun-like star (red), a lower main sequence star with a mass ~ 0.3Mq (purple), an upper main sequence star with a mass 
~ 7.9M0 (blue), a red super giant with a mass ~ 11.85M0 (orange), a white dwarf with a mass O.SMq (cyan), and a brown dwarf with 
a mass O.OSMq (brown). We plotted four curves for each star for comparison: the relativistic calculation for a co-rotating star (colored 
solid curve), the relativistic calculation for a non-rotating star (colored dashed curve), the Newtonian calculation for a co-rotating star 
(corresponding black solid curve), and the Newtonian calculation for a non-rotating star (corresponding black dashed curve). The vertical 
blue dashed lines represent the phase of inspiral due to gravitational radiation. 
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In order to understand the details of ALwind, a three- 
dimensional dynamical simulation will need to be employed. 
In this paper, for simplicity, we consider only the first two 
AL's, the mass accretion rate onto the hole will be the same 
as the stellar mass deposition rate. We then have: 



dL 
dt 



+ Tgr + m*jms = 0. 



(30) 



jma is the specific angular momentum at the ISCO, which 
equals in the Newtonian limit. 

For the same dwarf star with mass ~ O.SMq, we can 
calculate how its P and m change in the relativistic limit 
in comparison to the Newtonian treatment, as in Fig [T]and 
Fig. [51 With general relativity, the star does not recede from 
the hole as far as in the Newtonian limit, and the mass 
accretion rate has dropped to about half the Newtonian rate. 

We can apply the same method as in Section 2.3 to 
check if the stellar volume will e xceed the Roche volume 
in the evolution. Contradictory to teameurv et al. 1 l|l994l )'s 
result, the Roche evolution of a dwarf star is actually stable 
in the relativistic limit on dynamical time scales. 



4 CONSTRAINING SPINS FROM THE 
SIGNALS 

Returning to the observation of RE J 1034+396, the ~ Ihr 
periodic X-ray signal can be associated with a red giant 
or lower main-sequence star or brown dwarf orbiting and 
accreting onto a ~ 10® M© black hole or a ~ lO'^M© hole (the 
latter with a high spin). 16 cycles of signals are insufficient 
to be confident of a stable period. However, if more signals 
can be observed, information on spin could be extracted. 
In principle, the closer the star is to the hole, or the larger 
the spin of the hole, the stronger the emission would be. 
This could introduce an important selection effect as rapidly 
spinning holes would be more likely to be observed. 

Fig. llOl shows simulated periodic signals from a hot spot 
orbiting around the hole circularly using our numerical pack- 
age. The modulatio n of intensity com es from beaming effect 
- as we discussed in lDai et al.l ( 2010l ). the bolometric inten- 



sity varies as the fourth po wer of the Doppler-g ravitational 
redshift factor g = E/Eo (jCunninghaml [19751 ) . Here E is 
the energy of the photon in the observer's frame and Eq 
in the hot spot's inertial rest frame. The shape of the curves 
and the relative intensity are different for different spins and 
observer inclination angles. 

Also a major difference between producing a hot spot 
on the disc u sing the Roche accreti on model and other QPO 
models (e.g. [Wagoner et al.l (|200ll )) is that the periodicity 
of the modulation will follow that of the stellar orbit al- 
though there is likely to be phase noise. Since the stellar 
orbit evolves on an adiabatic time scale, the periodicity of 
the hot spot will also change very slowly but in a consistent 
way. For example, an accreting ~ O.3M0 dwarf star's orbital 
periodicity can change by as much as ~ Q.ls in one year, in 
which time we could observe ~ 8000 cycles of modulation. 
A P — P analysis similar to the double pulsar's could in 
principle be done to extract information of the stellar type 
and black hole mass and spin. 
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i=30 
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Figure 10. Simulated X-ray period emission from a hot spot or- 
biting a black iiole is shown for various black hole spin parameter 
a and observer inclination angle i. In each figure the hot spot 
orbits the hole in equatorial circular orbits with the same radius. 
The units are arbitrary but consistent. 



5 DISCUSSION AND FUTURE RESEARCH 
DIRECTIONS 

Observation of quasi-periodic X-ray emissions from stars or- 
biting black holes in Active Galactic Nuclei (AGN) is a po- 
tential probe of general relativity. There is an intriguing pos- 
sibility that the X-ray modulation of RE J1034-I-396 does 
come from an orbiting star filling its Roche lobe and be- 
ing swallowed b the black hole, i n a similar fashion to the 
" black widow" pulsar systems fe. g.[Thorsett fc Chakrabartvl 
J1999I)) like P S R195 7 +20 (e.g. iKulkarni fc Hester[ l|l988h : 
[Kulkarni et"al] (|l992l '): [Arons fc Tavanil (|l993h '). ItTs reason- 
able to search the existing AGN X-ray data base for periodic 
signals with relevant periods within range ~ 1 — 20hr. 

Although we are mostly interested in this paper in ob- 
servable tidally disrupted or accreting stars, it is quite rea- 
sonable to suppose that most dormant AGN black holes are 
orbited by a stable system of stars with a common orbital 
plane in much the same way as the sun is orbited by the 
planets. Sgr A*, which already has over 20 closely orbiting 
stars, might eventually evolve to such a configuration. Other 
orbiting stars can have also dynamical interactions with the 
accreting star and make its evolution more complicated. 

The probability of observing such phenomena relates 
to the stellar distribution in accretion discs and the mecha- 
nisms of bringing stars close to ISCO. When material accrete 
onto the black holes, stars can be borne in the accretion disc, 
which could go away eventually leaving only the stars. Such 
stars would orbit the hole in bound orbits in the equatorial 
plane with low eccentricity. On the other hand, if an accre- 
tion disc still exists when the star gets close to ISCO, the 
star will have interaction with the disc as well, or might be 
traveling in the disc already producing a gap and perturba- 
tion patterns on the disc. 

A series of other interesting future projects can also be 
done on this model. If the star is on a slightly eccentric or- 
bit or inclined orbit, the accretion rate could be larger than 
the case of a circular orbit, but perturbations will be added 
to the sinusoidal behavior of the fiux. A main-sequence star 
on its way to become a red-giant can also increase the mass 
accretion rate. A jet can also be formed due to the Roche ac- 
cretion - the radio component will also be very useful in pro- 
viding information on the black hole parameters. Between a 
total tidal disruption and a steady mass-transfer, interme- 
diate phases can also exist and we think much can be done 
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on that. In order to understand the stellar evolution and 
radiation mechanisms better, we also need to do a three- 
dimensional general relativistic dynamical simulation. The 
time it would take to run such evolutions on existing numer- 
ical codes would be too long. New numerical methods need 
to be designed. 

However, the most exciting and immediate prospect is 
to open up a new line of investigation using existing X-ray 
databases, to seek evidence for stars in bound orbits around 
massive black holes. Discovering such systems would elu- 
cidate the evolution of AGN and also open the possibility 
of a new measurement of black hole spin and ultimately of 
concluding that the Kerr metric correctly describes spinning 
black holes. 
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